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Abstract 

We study abelian dominance for confinement in terms of the local gluon 
properties in the maximally abelian (MA) gauge in a semi-analytical manner 
with the help of the lattice QCD. The global Weyl symmetry persistently 
remains as the relic of SU(A r c ) in the MA gauge, and provides the ambiguity 
on the electric and magnetic charges. We derive the criterion on the SU(iV c )- 
gauge invariance in terms of the residual symmetry in the abelian gauge. In 
the lattice QCD, we find microscopic abelian dominance on the link variable 
for the whole region of /3 in the MA gauge. The off-diagonal angle variable, 
which is not constrained by the MA-gauge fixing condition, tends to be 
random besides the residual gauge degrees of freedom. Within the random- 
variable approximation for the off-diagonal angle variable, we prove that 
off-diagonal gluon contribution to the Wilson loop obeys the perimeter law 
in the MA gauge, and show exact abelian dominance for the string tension, 
although small deviation is brought by the finite size effect of the Wilson 
loop in the actual lattice QCD simulation. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) is the fundamental theory of the strong interaction 
|]J2|]. Due to the asymptotic freedom, the gauge-coupling constant of QCD becomes small 
in the high-energy region and the perturbative QCD provides a direct and systematic 
description of the QCD system in terms of quarks and gluons. On the other hand, in 
the low-energy region, the strong gauge- coupling nature of QCD leads to nonperturbative 
features like color confinement, dynamical chiral-symmetry breaking and nontrivial 

topological effect by instantons and it is impossible to understand them directly 

from quarks and gluons in a perturbative manner. Instead of quarks and gluons, some 
collective or composite modes may be relevant degrees of freedom for the nonperturbative 
description in the infrared region of QCD. As for chiral dynamics, the pion and the sigma 
meson play the important role for the low-energy QCD, and they are included in the QCD 
effective theories like the (non-) linear sigma model, the chiral bag model PJTOfl and the 



Nambu-Jona-Lasinio model PJTTll, where these mesons are described as composite modes 
of quarks and anti-quarks. Here, the pion is considered to be the Nambu-Goldstone boson 



relating to spontaneous chiral-symmetry breaking and obeys the low-energy theorem [|12 
and the current algebra On the other hand, confinement is essentially described 
by dynamics of gluons rather than quarks. Hence, it is desired to extract the relevant 
collective mode from gluons for confinement phenomena. 

In 1970's, Nambu, 't Hooft and Mandelstam proposed an interesting idea that quark 
confinement can be interpreted using the dual version of the superconductivity [fT3|-|l5|l. 
In the ordinary superconductor, Cooper-pair condensation leads to the Meissner effect, 
and the magnetic flux is excluded or squeezed like a quasi-one-dimensional tube as the 
Abrikosov vortex Hl6 |, where the magnetic flux is quantized topologically. On the other 
hand, from the Regge trajectory [|17| of hadrons and the lattice QCD simulation p!8|JT9|] , 
the confinement force between the color-electric charge is characterized by the universal 
physical quantity of the string tension, and is brought by one- dimensional squeezing of 
the color-electric flux |?(J in the QCD vacuum. Hence, the QCD vacuum would be re- 
garded as the dual version of the superconductor based on above similarities on the low- 
dimensionalization of the quantized flux between charges. In this dual-superconductor 
picture for the QCD vacuum, the squeezing of the color-electric flux between quarks 
is realized by the dual Meissner effect as the result of condensation of color-magnetic 
monopoles, which is the dual version of electric-charge condensation. Monopole conden- 
sation and its relevant role for confinement were analytically pointed out by Seiberg and 
Witten very recently in the N = 2 supersymmetric version of QCD |2T| . 

However, there are two following large gaps between QCD and the dual superconductor 
picture. 
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1. This picture is based on the abelian gauge theory subject to the Maxwell-type 
equations, where electro-magnetic duality is manifest, while QCD is a nonabelian 
gauge theory. 

2. The dual-superconductor scenario requires condensation of (color-) magnetic 
monopoles as the key concept, while QCD does not have such a monopole as the 
elementary degrees of freedom. 

As the connection between QCD and the dual superconductor scenario, 't Hooft proposed 
the concept of the abelian gauge fixing EI|, the partial gauge fixing which only remains 
abelian gauge degrees of freedom in QCD. By definition, the abelian gauge fixing reduces 
QCD into an abelian gauge theory, where the off-diagonal element of the gluon field be- 
haves as a charged matter field and provides a color-electric current in terms of the resid- 
ual abelian gauge symmetry. As a remarkable fact in the abelian gauge, color-magnetic 
monopoles appear as topological objects corresponding to the nontrivial homotopy group 
n 2 (SU(A^ c )/U(l) JVc_1 ) = Z^" 1 . Thus, by the abelian gauge fixing, QCD is reduced into 
an abelian gauge theory including both the electric current j M and the magnetic current k^, 
which is expected to provide the theoretical basis of the monopole-condensation scheme 
for the confinement mechanism. 

As for the irrelevance of the off-diagonal gluon gluons, Ezawa and Iwazaki assumed 
abelian dominance that the only abelian gauge fields with monopoles would be essen- 
tial for the description of nonperturbative phenomena in the low-energy region of QCD, 
and showed a possibility of monopole condensation in the infrared scale by investigat- 
ing "energy-entropy balance" on the monopole current |23H26| in a similar way to the 
Kosterliz-Thouless transition in the 1+2 dimensional superconductivity f2"T| . Ezawa and 
Iwazaki formulated the dual London theory as an infrared effective theory of QCD, and 
later it is reformulated as the dual Ginzburg-Landau theory |2"5|-p3|]. 



Furthermore, abelian dominance [|35| -|3"9|j and monopole condensation |p3|-|42|j have been 
investigated using the lattice QCD simulation in the maximally abelian (MA) gauge 
4^-43|. The MA gauge is the abelian gauge where the diagonal component of the gluon 



is maximized by the gauge transformation. In the MA gauge, physical information of the 
gauge configuration is concentrated into the diagonal component as much as possible. The 
lattice QCD studies indicate abelian dominance that the string tension [j35| -|37| and the 
chiral condensate [|38|j39f| are almost described only by abelian variables in the MA gauge. 
Moreover, monopole dominance is also observed in the lattice QCD simulation in the MA 
gauge: only the monopole part in the abelian variable contributes to the nonperturbative 
QCD [p^,R3. Thus, the lattice QCD phenomenology suggests the dominant role of abelian 



variables including monopoles for the nonperturbative QCD in the MA gauge [ 44 ,45 



In this paper, we aim to understand the origin of abelian dominance for confinement 



3 



in terms of the local gluon properties in the MA gauge, and study the relation between 
macroscopic abelian dominance and microscopic abelian dominance in a semi-analytical 
manner with the help of the lattice QCD simulation. In Section 2, we study the residual 
symmetry and the gauge invariance condition for operators in the 't Hooft abelian gauge, 
paying the attention to the global Weyl symmetry. In Section 3, we investigate the MA 
gauge both in the lattice and in the continuum theories in terms of the gauge connection. 
In Section 4, we introduce the abelian projection rate as the overlapping factor between 
SU(2) and abelian link variables, and study microscopic abelian dominance in the MA 
gauge in the lattice formalism. In Section 5, we study the contribution of off-diagonal 
gluons to the Wilson loop in the MA gauge, and prove abelian dominance for the string 
tension in a semi-analytical manner. Section 6 is devoted to summary and concluding 
remarks. 



II. RESIDUAL SYMMETRY AND GAUGE INVARIANCE CONDITION IN 

THE ABELIAN GAUGE 

The dual superconductor picture for confinement phenomena is based on the abelian 



gauge theory including monopoles, and the 't Hooft abelian gauge fixing p2j is the key 
concept for the connection from QCD to such an abelian gauge theory. In this section, 
we investigate the abelian gauge fixing in QCD in terms of the residual gauge symmetry. 

The abelian gauge fixing, the partial gauge fixing which remains the abelian gauge 
symmetry, is realized by the diagonalization of a suitable SU(A c )-gauge dependent vari- 
able as $[A M (x)] G su(N c ) by the SU(iV c ) gauge transformation. In the abelian gauge, 
^[A^x)] plays the similar role of the Higgs field, and can be regarded as the composite 
Higgs field. 

For an hermite operator $[A M (x)] which obeys the adjoint transformation, is 
transformed as 

$(z) = $ a T a = n(x)<S>(x)&(x) = H ■ $ diag (x) 

= diag(A 1 (x),---,A 7Vc (x)), (1) 

using a suitable gauge function Q(x) = exp{i£ a (x)T a } G SU(iV c ). Here, each diagonal 
component X 1 (i=l,- ■ -,N C ) is to be real for the hermite operator $[A M (x)]. In the abelian 
gauge, the SU(A C ) gauge symmetry is reduced into the U(l) 7Vc ~ 1 gauge symmetry corre- 
sponding to the gauge-fixing ambiguity. The operator is diagonalized to H ■ <& d ia g {x) 
also by the gauge function fl?(x) = u(x)Q(x) with uj(x) = exp(—iH ■ (p{x)) G U(l) JVc_1 , 

-> n"(x)$(x)Q^(x) = uj(x)H ■ $ diag (x)J(x) = H ■ $ diag (x), (2) 
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and therefore U ( 1 ) ^ — 1 abelian gauge symmetry remains in the abelian gauge. 

In the abelian gauge, there also remains the global Weyl symmetry as a "relic" of 
the nonabelian theory [|46 , [47| . Here, the Weyl symmetry corresponds to the subgroup of 



SU(iV c ) relating to the permutation of the basis in the fundamental representation. Then, 
the Weyl group is expressed as the permutation group Pjv c including Ar c C 2 elements. 
For simplicity, let us consider the N c = 2 case. For SU(2) QCD, the Weyl symmetry 
corresponds to the interchange of the SU(2)-quark color, |+) = (J) and |— ) = (5), in the 
fundamental representation. The global Weyl transformation is expressed by the global 
gauge function, 

W = ^cosa+f sina^r = CQg Q + r2sina } 

/ o e~ ia \ 

= <( e *a o j eP 2 GSU(2) (3) 

with an arbitrary constant a G R. By the global Weyl transformation W, the SU(2)-quark 
color is interchanged as = ie ia \— ) and W\— ) = ie~ ia \+) except for the global phase 

factor. This global Weyl symmetry remains in the abelian gauge, because the operator 
is also diagonalized by using Q w (x) = Wfl(x), 

•(,) - n-www = <"W*)i«* = -w*)?. (4) 

Here, the sign of &diag(x), or the order of the diagonal component X l (x), is globally 
changed by the Weyl transformation. It is noted that the sign of the U(l) 3 gauge field 
= is globally changed under the Weyl transformation, 

A ll -+A™ = WA^Wt = -4| = -A,. (5) 

Therefore, all the sign of the abelian field strength, electric and magnetic charges are also 
globally changed: 

JT = P Z£ i -rVK xxr <r rxr\ -r" 

J flU — [W 2 /ij/ YY J ^ V VY J ^is, 

kjj, = d J~afj, * k„ = k/j,. (6) 



In the abelian gauge, the absolute signs of the electric and the magnetic charges are 
settled, only when the Weyl symmetry is fixed by the additional condition. When $L4 M (;r)] 
obeys the adjoint-type gauge transformation like the nonabelian Higgs field, the global 
Weyl symmetry can be easily fixed by imposing the additional gauge-fixing condition as 
&diag(x) > for SU(2), or the ordering condition of the diagonal components A* in H-^diag 
as A 1 > .. > X Nc for the SU(iV c ) case. As for the appearance of monopoles in the abelian 
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gauge, the global Weyl symmetry Pjv c is not relevant, because the nontriviality of the 
homotopy group is not affected by the global Weyl symmetry. However, the definition of 
the magnetic monopole charge, which is expressed by the nontrivial dual root of SU(iV c ) dual 
23| , is globally changed by the Weyl transformation. 

Now, we consider the abelian gauge fixing in terms of the coset space of the fixed gauge 
symmetry. The abelian gauge fixing is a sort of the partial gauge fixing which reduces 
the gauge group G =SU(iV c )i oca i of the system into its subgroup H =U(1) j^cai 1 ( x Pf/° bal ) 
including the maximally torus subgroup of G. In other words, the abelian gauge fixing 
freezes the gauge symmetry relating to the coset space G/H, and hence the representative 
gauge function Q which brings the abelian gauge belongs the coset space G/H: Q G G/H. 
In fact, Q G G/H is uniquely determined without the ambiguity on the residual symmetry 
H, if the additional condition on H is imposed for Q. 

However, such a partial gauge fixing makes the total gauge invariance unclear. Here, 
let us consider the SU(iV c ) gauge- invariance condition on the operator defined in the 



abelian gauge f4q| . To begin with, we investigate the gauge-transformation property of 
the gauge function Q G G/H which brings the abelian gauge (See. Fig.l). For simplicity, 
the operator $ to be diagonalized is assumed to obey the adjoint gauge transformation 
as $ — > $f = g$gt with y g G G. After the gauge transformation by y g G G, Vt 9 G G/H 
is defined so as to diagonalize $ 9 as Q 9 § 9 (Q 9 Y = &diag, and hence the gauge function 
fl 9 G G/H which realizes the abelian gauge is transformed as 

Q -> Q 9 = h[g]VLg ] (7) 

under arbitrary SU(iV c ) gauge transformation by g G G. Here, h[g] G H is chosen so as 
to make Q 9 belong G/H, i.e., Q 9 G G/H. (If the additional condition on H is imposed to 
specify Q G G/H, Qg^ does not satisfy it in general.) This is similar to the argument on the 
hidden local symmetry [|l^] in the nonlinear representation. In general, the gauge function 
Q G G/H is transformed nonlinearly by the gauge function g due to h[g] G H. Thus, the 
gauge-transformation property on the gauge function Q G G/H becomes nontrivial in the 
partial gauge fixing. 

Owing to the nontrivial transformation (|7]) of Q &G/H, any operator O n defined in 
the abelian gauge is found to be transformed as O n — > (O^)^ by the SU(iV c ) gauge 
transformation of v g G G. We demonstrate this for the gluon field = Q(A^ + ^8^)0) 
in the abelian gauge. By the gauge transformation of ^ g G G, is transformed as 



% 

—i 

e 

Here, we have used 

-dp)g\ = {g 2 gi){A tl + -, 

e e 



A° - (Af.r = A™ 9 = A^ n = (Aj) fc W = %]« + -d,W[g\. 



KT 2 = 92{Af l 1 + \)gl = (g 2 g 1 )(A tl + -d,){g 2 g^ = (A^) 9291 (9) 
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for the successive gauge transformation by gi, g 2 G SU(iV c ). Similarly, the operator O n 
defined in the abelian gauge is transformed by v g G G as 

n _^ (Q B yv> = Q9 9Q^ = h[g]VLg ] ■ gOg ] ■ g^h)[g] 

= h\g]DOQW\g] = h[g}O n rf[g] = (O a ) h ^\ (10) 

as shown in Fig.l. Here, O is assumed to obey the adjoint transformation as O 9 = gOg^ 
for simplicity. 

Thus, arbitrary SU(iV c ) gauge transformation by g G G is mapped into the partial 
gauge transformation by h[g] G H for the operator O n defined in the abelian gauge, and 
O n transforms nonlinearly as O n — > (O^)^ 9 ' by the SU(iV c ) gauge transformation g. If 
the operator O n is if-invariant, one gets (O a ) h ^ = O n for any h[g) G H, and hence 
O is also G-invariant or total SU(iV c ) gauge invariant, because O is transformed into 
(0°)%] = O n by y g G G. Thus, we find a useful criterion on the SU(iV c ) gauge mvanance 
of the operator defined in the abelian gauge j|6| : If the operator O defined in the abelian 
gauge is if-invariant, O n is also invariant under the whole gauge transformation of G. 

Here, let us consider the application of this criterion to the effective theory of QCD 
in the abelian gauge, the dual Ginzburg-Landau (DGL) theory p8| , |29|| . In the DGL 
theory, the local U(l) 7Vc_1 and the global Weyl symmetries remain, and the dual gauge 
field and the monopole field Xa [ck=1, ■ ■ •, ^N C (N C — 1)] are the relevant modes for 
infrared physics. Although is invariant under the local transformation of U(l) iVc 1 C 
SU(iV c ), B^ = B^ ■ H is variant under the global Weyl transformation, and therefore _B M 
is SU(iV c )-gauge dependent object and does not appear in the real world alone. As for 
the monopole field, there exists one Weyl-invariant combination of the monopole field 
fluctuation, x = J2 a Xa p2| , which is also U(l) 7Vc_1 -invariant. Therefore, the monopole 
fluctuation x is completely residual-gauge invariant in the abelian gauge, so that x is 
SU(A c )-gauge invariant and is expected to appear as a scalar glueball with J PC = ++ , 
like the Higgs particle in the standard model. 



III. MAXIMALLY ABELIAN (MA) GAUGE IN THE CONNECTION 

FORMALISM 

The abelian gauge has some arbitrariness corresponding to the choice of the operator $ 
to be diagonalized. As the typical abelian gauge, the maximally abelian (MA) gauge, the 
Polyakov gauge and the F12 gauge have been tested on the dual superconductor scenario 
for the nonperturbative QCD. Recent lattice QCD studies show that infrared phenomena 
such as confinement properties and chiral symmetry breaking are almost reproduced in 



the MA gauge |35|-^T|. In the SU(2) lattice formalism, the MA gauge is defined so as to 
maximize 
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Rma[U,] = £tr{^( S )T 3 ^( S )T 3 } 

S,fJ, 

= 2 EK(*) 2 + Ul{sf - U^sf - Ul{sf} 

S,fL 

= 2^[i-2{(/;( s f + f/;( s ) 2 }] (ii) 



by the SU(2) gauge transformation. Here, the link variable U^(s) = U®(s) + ir a U®(s) G 
SU(2) with E/J(s), G R relates to the (continuum) gluon field A M = A a ^T a G 

sit (2) as t^(s) = e Mej4 ''^ s \ where e denotes the QCD gauge coupling and a the lattice 
spacing. In the MA gauge, the absolute value of off-diagonal components, U^(s) and 
U^(s), are forced to be small. In the continuum limit a — > 0, the link variable reads 
U^(s) = e weA M( s ) — i _)_ iaeA M (s) + 0(a 2 ), and hence the MA gauge is found to minimize 
the functional, 

R ch [A,} ee ie 2 J d'x{Al(xf + A%xf} = e 2 J d*xA+(x)A;(x), (12) 

with Ap(x) = A={A^(x) ± iy4 2 (a;)}. Thus, in the MA gauge, the off-diagonal gluon 
component is globally forced to be small by the gauge transformation, and hence the 
QCD system is expected to be describable only by its diagonal part approximately. 

The MA gauge is a sort of the abelian gauge which diagonalizes the hermite operator 

^K(s)] = J2U±^s)r 3 uUs). (13) 

Here, we use the convenient notation U^^(s) = Uj t (s — ft) in this paper. In the continuum 
limit, the condition of the MA gauge becomes £(i9 M ^ e Af l )A^ = 0- This condition can 

be regarded as the maximal decoupling condition between the abelian gauge sector and 
the charged gluon sector. 

In the MA gauge, $(s) is diagonalized as $ma(s) = ^dia g {s)r 3 with &di ag (s) G R-, 
and there remain the local U(l)3 symmetry and the global Weyl symmetry f4"6fl . As a 
remarkable fact, $(s) does not obey the adjoint transformation in the MA gauge, and the 
sign of &diag(s) is not changed by the Weyl transformation by W in Eq.@, 

*MA(s) = ®diag(s)T 3 

= -E WU ± ,(s)r 3 Ul,(s)W^ = -W^ diag {s)r 3 W^ = <5> dmg (s)r 3 . (14) 

Thus, the Weyl symmetry is not fixed in the MA gauge by the simple ordering condition 
as Qdiag > 0, unlike the adjoint case. We find the gauge invariance condition on the 
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operator O n defined in the MA gauge: if O n is invariant both under the local U(l) Ar<; ~ 1 
gauge transformation and the global Weyl transformation, O n is also invariant under the 
SU(iV c ) gauge transformation. 

In the continuum SU(iV c ) QCD, it is more fundamental and convenient to define the 
MA gauge fixing by way of the SU(A" c )-covariant derivative operator = <9 M + ieA^, 
where is the derivative operator satisfying [<9 M , f(x)] = d^f(x). The MA gauge is 
defined so as to make SU(iV c ) -gauge connection = <9 M + ieA^T a close to U(l) Arc ~ 1 - 
gauge connection D^ hcl = <9 M + ieA^ ■ H by minimizing 

R ch = J d A x tr[£„, H}=e 2 J d A x tr[A M , H]^, H] 

, p 2 . N C (N C -1) 

= e 2 jd*xY.K*K*-MElE p ) = e -jd*x Yl K\ 2 > ( 15 ) 

a,/3 a=l 

which expresses the total amount of the off-diagonal gluon component. Here, we have used 

N C (N C -1) 

the Cartan decomposition, = A a ^T a = A^-H+ ^ A°E a ; H = (T 3 , T 8 , • • • , T N 2_ Y ) 

a=l 

is the Cartan subalgebra, and E a (a = 1, 2, • • • , N 2 — N c ) denotes the raising or lowering 
operator. In this definition with D^, the gauge transformation property of R c h becomes 
quite clear, because the SU(iV c ) covariant derivative obeys the simple adjoint gauge 
transformation, — > VLD^Vtf, with the SU(iV c ) gauge function VL G SU(iV c ). By the 
SU(A" C ) gauge transformation, R c ^ is transformed as 

R ch -> R^ = J d A x tr (pDfl, H^[nD^, H]) 

= J d A x tr ([£>„, fi t #fi] t [ J D M , Q^HQ]) , (16) 

and hence the residual symmetry corresponding to the invariance of R^ is found to be 
U(l)^ 1 x P g ^ bal cSU(Ay iocal , where P g ^ bal denotes the global Weyl group relating to 
the permutation of the N c basis in the fundamental representation. In fact, one finds 
lu^Hlu = H for uj = e -^ x > S G U(l)^~ 1 1 , and the global Weyl transformation by W G 
-fgiobai om y exchanges the permutation of the nontrivial root dj and never changes P c h- In 
the MA gauge, arbitrary gauge transformation by v f2 G G is to increase R ch as R^ h > R c h- 
Considering arbitrary infinitesimal gauge transformation Q = e l£ ~ 1+ie with v e Esu(N c ), 

— * — * — * 

one finds Q*HQ ~ H + e] and 

P c n h ~ R ch + 2if d 4 xtr ([£)„, [H, e]} [P> M , H]) 

= R ch + 2iJ d A xti (e[H, [£>„, [£>„, #]]]) . (17) 
In the MA gauge, the extremum condition of R^ h on V £ Esu(N c ) provides 

[H,[D„[D„H]]] = 0, (18) 
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which leads to ± eAfyAf^ = for the N c =2 case. Thus, the operator $ to be 

diagonalized in the MA gauge is found to be 

$IA,] = [D„[D„H}} (19) 

in the continuum theory. 



IV. MICROSCOPIC ABELIAN DOMINANCE IN THE MA GAUGE 



In the abelian gauge, the diagonal and the off-diagonal gluons play different roles in 
terms of the residual abelian gauge symmetry: the diagonal gluon behaves as the abelian 



gauge field, while off-diagonal gluons behave as charged matter fields ||40|1 . Under the 
U(l) 3 gauge transformation by uj = exp(— i<£>y) G U(l) 3 , one finds 

4 - W = A l + ( 20 ) 

K - (^) ± = A%e^ (21) 

with A^ = -^(^ ±iA^). The abelian projection is simply defined as the replacement of 
the gluon field A^ = A aJ ^ G su(2) by the diagonal part = A 3 ^ G u(l) 3 C su(2). 

We call "abelian dominance for an operator 0[yl M ]", when the expectation value (O) 
is almost unchanged by the abelian projection A^ — > as (OL4.J) ~ (0[A^])a.g., when 
( )a.g. denotes the expectation value in the abelian gauge. Ordinary abelian dominance 
is observed for the long-distance physics in the MA gauge, and this would be physically 
interpreted as the effective-mass generation of the off-diagonal gluon induced by the MA 
gauge fixing [Ti8|Ji9[ . 



In the lattice formalism, the SU(2) link- variable U^s) can be factorized as 



U^s) = M^s)u^s) G G 



M,(s) = exp (i{Ti0j( S ) + r 2 el(s)}) G G/H, 
u^s) = exp (ir 3 6l(s)) G H (22) 

with respect to the Cartan decomposition of G = G/H x H into G/H =SU(2)/U(1) 3 
and H =U(1) 3 . Here, the abelian link variable, 








^ (s) = n E U(l) 3 C SU(2), (23) 



plays the similar role as the SU(2)-link variable U^(s) G SU(2) in terms of the residual 
U(l) 3 gauge symmetry in the abelian gauge, and 0^(s) G (— 7r,7r] corresponds to the 
diagonal component of the gluon in the continuum limit. On the other hand, the off- 
diagonal factor M M (s) gSU(2)/U(1) 3 is expressed as 
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/ cos^(s) -sin^(s)e- ix ^ s ) 
~ \fa.n 6 li (s)jx»W cosfl M (s) 
-|c,( S )| 2 -<£(*) \ 



(24) 



00 v /rr M^ 



with On(s) = mod^y (6 1 *) 2 + (# 2 ) 2 G [0, |] and X/i( s ) G ( — t 1 ") 71 "]- Near the continuum 
limit, the off-diagonal elements of M M (s) correspond to the off-diagonal gluon components. 
Under the residual U(l)3 gauge transformation by u(s) = e~ vp ( a )~* G U(l) 3 , u^(s) and 
M M (s) are transformed as 

M s ) ^«£( s ) =u{s)u li (s)ui(s + fl) EH (25) 
M^s) - M;( S ) = W ( S )M M ( S )o;t( s ) g (26) 

so as to keep M£(s) belong G/H. Accordingly, 6^(s) and c M (s) G C are transformed as 

*J00 - 0?OO = mod 27r [^( S ) + + //) - v7( a )}/2] (27) 
c«W - = c„( a )e*« (28) 



Thus, on the residual U(l) 3 gauge symmetry, w M (s) behaves as the U(l) 3 lattice gauge 
field, and 0^(s) behaves as the U(l) 3 gauge field in the continuum limit. On the other 
hand, M jU (s) and c M (s) behave as the charged matter field in terms of the residual U(l) 3 
gauge symmetry, which is similar to the charged weak boson in the standard model. 

In this parameterization (|22"D, there are two U(l)-structures embedded in SU(2) cor- 
responding to e ld v and e l * M . To clarify this structure, we reparametrize the SU(2) link 



variable as 



. N / cos9 u e l9 ^ — sin^e \ , . 

M ; \sm9„e 1 ^ cos^e-< / V ; 



or equivalently 



u l = cos °n cos e l> u l = sin Op sin x M , 

U* = cos 6> M sin 0j, f/ 2 = sin 6> M cos x M , (30) 

with Xfj, = X/x + The range of the angle variable can be redefined as < # M < | and 
— 7r < O^iXn — 71 "■ Here, (U^U^U^U 3 ) forms an element of the 3-dimensional hyper- 
sphere S 3 ~SU(2), because of (f/°) 2 + (^) 2 + (Uff + (U 3 ) 2 = 1. For a fixed 0„, both 
(U®, U 3 ) and ([/* £/ 2 ) form the two 5 1 ~U(l) subgroups embedded in S 3 in a symmetric 
manner. From the parametrization in Eq. (|29"D , the SU(2) measure can be easily found as 
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/ dU, = f dUXUffi 5 ( £ ( Up 2 - 1) 

J J a=0 

= A [* ^sin^cos^ f dx„ T ^J- (31) 

In the lattice formalism, the abelian projection is defined by replacing the SU(2) link 
variable U^(s) G SU(2) by the abelian link variable u^(s) G U(l) 3 . 

In the MA gauge, the off-diagonal gluon component is strongly suppressed, and the 
SU(2) link variable is expected to be U(l)3-like as U^s) ~ u^(s) in the relevant gauge 
configuration. In the quantitative argument, this can be expected as (f/ At (s)M| t (s))MA — 1, 
where ( )ma denotes the expectation value in the MA gauge. In order to estimate the 
difference between t/ M (s) and u^(s), we introduce the "abelian projection rate" i^Abei 
50| , |5T|] , which is defined as the overlapping factor as 

RAbei(s,/i) = ^Re tr{Up(s)ul(s)} = ^Re trM„(s) = cos^(s). (32) 

This definition of i?Abei is inspired by the ordinary "distance" between two matrices 
A,B G GL(iV,C) defined as d 2 (A, B) = |tr{(A - B)\A - B)} [||, which leads to 
d 2 {A,B) = 2 - Re tr(ABt) for A, B gSU(2). The similarity between U^s) and u^s) 
can be quantitatively measured in terms of the "distance" between them. For instance, 
if cos 6^(5) = 1, the SU(2) link variable becomes completely abelian as 



, e-^ J 1 

while, if cos 6^(5) = 0, it becomes completely off-diagonal as 








We show in Fig. 2(a) the probability distribution P(i?Abei) of the abelian projection 
rate i?Abei(s,/x) = cos 9^(s) in the MA gauge. Here, (RAbei) 0=0 m the strong coupling 
limit (P = 0) |50| , |5T| is analytically calculable as 



/p \ _ /msfl / e « _ I dU^{s) cos e^s) 



JdU,(s) 
J f dO^s) sinO^s) cos 2 O^s) _ 2 
dO^s) sinO^s) cos O^s) 3 



(33) 



using Eq. (|3lD . In the MA gauge, i?Abei approaches to unity as shown in Fig.2(a). The 
off-diagonal component of the SU(2) link variable is forced to be reduced. As a typical 
example, one obtains (-RAbei)MA — 0.926 on 16 4 lattice with j3 = 2.4. We show also 
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the abelian projection rate (-RAbei)MA as the function of (3 in Fig.2(b). For larger /3, 
(cos ^(s)}ma becomes slightly larger. Without gauge fixing, the average (i?Abei) is found 
to be about | without dependence on (3. In the continuum limit in the MA gauge, 
U*(s) and U 2 (s) become at most 0(a), and therefore (i?Abei)MA approaches to unity 
as (i?Abei)MA = 1 + 0(a 2 ) due to the trivial dominance of U°(s), which differs from 
abelian dominance in the physical sense. The remarkable feature of the MA gauge is the 
high abelian projection rate as (-RAbei)iviA — 1 in the whole region of (3. In fact, we find 
(i?Abci)MA — 0.88 even for the strong coupling limit (3 = 0, where the original link variable 
Ufj, is completely random. Thus, abelian dominance for the link variable is observed 
at any scale in the MA gauge. 

To understand the origin of the high abelian projection rate as (i?Abei)MA — 1, we 
estimate the lower bound of (-RAbei)MA m the MA gauge using the statistical consideration. 
The MA gauge maximizes 

iWtg = J2^{ U ^Ul(s)r 3 } = tr(T 3 £0 M ( S )), (34) 

where M (s) = U^s^UUs) is an su(2) element satisfying 0^ = 1. Denoting M (s) = 
0^(s)r a , we parameterize the 3-dimensional unit vectors M = (0^,0^, 0^) G S 2 (// = 
1,2,3,4) as M = (sin 29^ cosx^, sin 29 ^ sin x M , cos 29^) using Eqs.(f22j) and The MA 
gauge maximizes the third component 0^ using the gauge transformation. Under the local 
gauge transformation by V(s) = 1 + {V(s ) — 1}5 SS0 G SU(2), (t)^(s ) is transformed as 
the unitary transformation, 

Mso) - ^(so) = V(s )Mso)V~\s ), (35) 

which leads to a simple rotation of the unit vectors M . In the MA gauge, ^0^ is "po- 

S,fl 

larized" along the positive third direction. On the 4-dimension lattice with N sites, 4A" 
unit vectors (f>n(s) are maximally polarized by N gauge functions V(s) in the MA gauge. 
Then, (i?Abei)MA is expressed as the maximal "polarization rate" of 4A" unit vectors M by 
suitable N gauge functions V(s). On the average, this estimation of (-RAbei)MA is approx- 
imately given by the estimation of the maximal polarization rate of 4 unit vectors M by 
a suitable rotation with V G SU(2). The lower bound of (-RAbei)MA is obtained from the 
strong-coupling system with (3 = 0, where link variables U^s) are completely random. 
Accordingly, M can be regarded as random unit vectors on S 2 . The maximal "polariza- 

4 4 

tion" of 4 unit vectors M is realized by the rotation which moves =X]0/x/lX/ < <M ^° ^ ne 

unit vector cf) R = (0,0, 1) in third direction. Here, cos 29 R after the rotation is identical 
to the inner product between M and 0, because of • M = (fi R ■ 4> R = (0^) 3 = cos 29 R . 
Then, we estimate (i^Abei)MA = (cos# m )ma at (3 = as 
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(cos0X= A ° ^{{[ / dU,} (- £ cos^ 

= ill - / 2 d6^ cos 6^ sin 8^ \ dx»} 7 51 cos { 9 cos _1 (0- M )} . (36) 

Using this estimation (|36|) , we obtain (-RAbei)MA — 0.844, which is close to the lattice 
result (-RAbei) — 0.88 in the strong coupling limit (/3 = 0). Such a high abelian rate 
(R Abel) ma in the MA gauge would provide a microscopic basis of abelian dominance for 
the infrared physics. 

V. SEMI- ANALYTICAL PROOF OF ABELIAN DOMINANCE FOR 

CONFINEMENT FORCE 

Abelian dominance and monopole dominance for the nonperturbative phenomena are 
numerically observed in the MA gauge in the lattice QCD simulations [^,^,0-49]. As 



for confinement, monopole dominance, which means the dominant role of the magnetic 
current k^ than that of the electric current j M , seems trivial if abelian dominance holds, 
because does not provide the electric confinement in 1+3 dimension. Then, as for 
confinement, abelian dominance, which means the dominant role of the diagonal element 
than that of the off-diagonal element, is the nontrivial interesting phenomenon observed 
in the MA gauge in the lattice QCD. In this section, we study the origin of abelian 
dominance on the string tension as the confinement force in a semi-analytical manner, 
considering the relation with microscopic abelian dominance on the link variable [[51 . 



In the MA gauge, the diagonal element cos^(s) in M^s) is maximized by the gauge 
transformation as large as possible. For instance, the abelian projection rate is almost 
unity as i?Abei = (cos^(s))ma — 0.93 at f3 = 2.4. Then, the off-diagonal element 
e *XA«0) sin^ At (s) is forced to take a small value in the MA gauge due to the factor sin^(s), 
and therefore the approximate treatment on the off-diagonal element would be allowed in 
the MA gauge. Moreover, the angle variable Xn( s ) i> s n °t constrained by the MA gauge- 
fixing condition at all, and tends to take a random value besides the residual U(l) 3 gauge 
degrees of freedom. Hence, Xni s ) can be regarded as a random angle variable on the 
treatment of M^(s) in the MA gauge in a good approximation. 

Let us consider the Wilson loop (W c [U^(s)]) = (trjjt^(s)) = (tr JJ{M jU (s)M M (s)}) 

c c 

in the MA gauge. In calculating (WciU^s)}), the expectation value of e %x ^ in M M (s) 
vanishes as 



2tt 



(e*M«) ~ / d X ^s) exp{^(s)} = 0, (37) 



14 



when Xni s ) behaves as a random angle variable. Then, within the random- variable ap- 
proximation for the off- diagonal factor M (U (s) appearing in (yV c\U ^(s)}) is simply 
reduced as a c-number factor, M^s) — > cos 9^(s) 1, and therefore the SU(2) link variable 
Un(s) in the Wilson loop {Wc\U ^(s)]) is simplified as a diagonal matrix, 

U^s) = M^u^s) -> cos^(s)^(s). (38) 

Then, for the I x J rectangular C, the Wilson loop Wc[£^(s)] in the MA gauge is 
approximated as 

(WcWM) = (trllM*)) - (II 008 M*<) ^11 ««( fl i)>MA 

i=l i=l j=l 

L 

~ (exp{^m(cos0 w (sj))}) M A (W / cK( s )])ma 
~exp{L(ln(cos^(s))) M A} (W c K{s)}) MA , (39) 

L 

where L = 2(1+ J) denotes the perimeter length and Wcftt^s)] = tr JJ the abelian 

i=i 

L 

Wilson loop. Here, we have replaced ^ln{cos(0 Mi (si)} by its average L(ln{cos6' A1 (s)})MA 

i=l 

in a statistical sense, and such a statistical treatment becomes more accurate for larger 
/, J and becomes exact for infinite /, J. 

In this way, we derive a simple estimation as 

ee (W c [[/ m ( S )])/(W c [^( S )]}ma ^ exp{L(ln(cos# M ( S ))>MA} (40) 

for the contribution of the off-diagonal gluon element to the Wilson loop. From this 
analysis, the contribution of off-diagonal gluons to the Wilson loop is expected to obey 
the perimeter law in the MA gauge for large loops, where the statistical treatment would 
be accurate. 

Now, we study the behavior of the off-diagonal contribution 
Wq S ee (W C [U M (s)]) / (W c [u^(s)}) ma in the MA gauge using the lattice QCD, considering 
the theoretical estimation Eq.(^). As shown in Fig.3, we find that seems to obey 
the perimeter law for the Wilson loop with /, J > 2 in the MA gauge in the lattice QCD 
simulation with f3 = 2.4 and 16 4 . We find also that the behavior on Wq S as the function of 
L is well reproduced by the above analytical estimation with the microscopic information 
on the diagonal factor cos^(s) as (ln{cos6' /J (s)})MA — —0.082 for (3 = 2.4. Thus, the 
off-diagonal contribution Wq to the Wilson loop obeys the perimeter law in the MA 
gauge, and therefore the abelian Wilson loop {Wc[u^(s)}) ma should obey the area law as 
well as the SU(2) Wilson loop Wc[t^(s)]. From Eq. (flop , the off-diagonal contribution to 
the string tension vanishes as 
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A(T = CT S U(2) - OAbel (41) 

= ~ ijwoo W ln ( W RMU^s)}) + R jxni oo \n(W RxT [u^s)}) M A 

R-\-T 

~ -2(ln{cos^(s)}) M A Jim -==- = 0. 

Thus, abelian dominance for the string tension, ctsu(2) = 0"Abei> can be proved in the MA 
gauge by replacing the off-diagonal angle variable x^{ s ) as a random variable. 

The analytical relation in Eq.(M) indicates also that the finite size effect on R and 
T in the Wilson loop leads to the deviation between the SU(2) string tension ctsu(2) and 
the abelian string tension CTAbci as ctsu(2) > 0"Abei in the actual lattice QCD simulations. 
Here, we consider this deviation Act = ctsu(2) — o"Abei in some detail. Similar to the SU(2) 
inter-quark potential Vsu(2)( r ) from ( Wsu(2) ) = ( )> we define the abelian inter- 

quark potential VAbei(^) and the off-diagonal contribution V s(r) of the potential from 
( W A be\ ) = ( W[u^(s)] ) and W oS , respectively, 

Vsm (r) = -^-ln{W sm (RxT)), 

V Abel (r) = ~ In ( W Ahel (R x T) ), 
1 a 

V oS (r) = -jrJnWMR * T) = --In ( W ^ {R - T) ) 

= V su(2) (r)-y Abel (r), (42) 

wehre r = Ra denotes the inter-quark distance in the physical unit. We show in Fig. 4 
^su(2)( r ), VAbei( r ) and V Q g(r) extracted from the Wilson loop with T = 7 in the lattice 
QCD simulation with (3 = 2.4 and 16 4 . As shown in Fig. 4, the lattice result for V g(r) 
seems to be reproduced by the theoretical estimation obtained from Eq.(|40|), 

Kfr(r) = V sm) (r) - V Ahel {r) ^ _ 2 (^+ T ) ( ln(cos ^( g ))) MA (43) 

using the microscopic information of (1ii(cos^ /1 (s)))ma — —0.082 at (3 = 2.4. From the 
slope of V g(r) in Eq . (|43|) , we can estimate Act = ctsu(2) — cr Abel in the physical unit as 

Act = ct su(2) - CT Abcl ~ -2( In ( cos 6> M (s) ) ) MA 



To? 



= -(ln(l-sin 2 ^( S ))) MA ^. (44) 

In the MA gauge, sin 2 # M (s) takes a small value and can be treated in a perturbation 
manner so that one finds 

Act ~ (sin 2 M ( a )>MA-^ = <(^(*)) 2 + K 2 {s)Y) UA ^- 2 . (45) 
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Near the continuum limit a ~ 0, we find U£ ~ aeA^/2 (a=l,2,3) from = e iaeA ^ T " ' I 2 , 
and then we derive the relation between Act and the off-diagonal gluon in the MA gauge 
as 

Aa ~ ^({eAlf + (e4) 2 >MA = ^(Mj) 2 + M^W, (46) 

where t = Ta is the temporal length of the Wilson loop in the physical unit. In Eq . (f46|) , 
{(eA 1 ^) 2 + (eA 2 t ) 2 )yix is the off-diagonal gluon-field fluctuation, and is strongly suppressed 
in the MA gauge by its definition. It would be interesting to note that microscopic 
abelian dominance or the suppression of off-diagonal gluons in the MA gauge is directly 
connected to reduction of the deviation Act in Eq.(fK]). Since ((eA 1 ^) 2 + (eA^) 2 ) M A is a 
local continuum quantity, it is to be independent on both a and t. Hence, the deviation 
Aa between the SU(2) string tension <7gu(2) and the abelian string tension CTAbei can be 
removed by taking the large Wilson loop as t —>■ oo or the small mesh as a —>■ with fixed 
t. 

Finally in this section, we study the origin of abelian dominance in the MA gauge in 
terms of the properties of the off-diagonal element 

c,(s)=e ix ^ s) sm9,(s) (47) 

of M M (s) in the link variable U^s), considering the validity of the random- variable ap- 
proximation for Xfj,{ s ) i n the MA gauge. In the above treatment, the contribution of the 
off-diagonal element in the link variable U^s) is completely dropped off, and its effect 
indirectly remains as the appearance of the c-number factor cos6 l At (s) in the link variable. 
Such a reduction of the contribution of the off- diagonal elements is brought by the two 
relevant features on the two local variables, O^s) and x M (s), in the MA gauge. One is the 
microscopic abelian dominance as (cos# m (s))ma — 1 in the MA gauge, and the other is 
the randomness of the off-diagonal variable Xm( s )- 

1. In the MA gauge, the microscopic abelian dominance holds as (cos^ a1 (s))ma — 1, 
and the absolute value of the off-diagonal element |c M (s)| = | sin ^(s) | is strongly 
reduced. Such a tendency becomes more significant as (3 increases. 

2. The off-diagonal angle variable x M (s) is not constrained by the MA gauge-fixing 
condition at all, and tends to be a random variable. In fact, x^i s ) is affected 
only by the action factor e _/35QCD in the QCD generating functional, but the effect 
of the action to Xm( s ) is quite weaken due to the small factor sin 6^(3) in the MA 
gauge. The randomness of x^{ s ) tends to vanish the contribution of the off-diagonal 
elements. 
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Here, the randomness of the off-diagonal angle- variable x^i s ) is closely related to the 
microscopic abelian dominance. In fact, the randomness of Xv( s ) is controlled only by 
the action factor e _/3SQCD in the QCD generating functional, however the effect of the 
action to i s quite weaken due to the small factor sin^ At (s) in the MA gauge, because 

Xfi(s) always accompanies sin^(s) in the link variable U^(s). Near the strong-coupling 
limit (3 ~ 0, the action factor e _/35QCD brings almost no constraint on X/t( s ) m the MA 
gauge. The independence of Xni s ) from the action factor is enhanced by the small factor 
sin6 l M (s) accompanying x^{ s )- Hence, behaves as a random angle-variable almost 

exactly, and the contribution of the off-diagonal element is expected to disappear in the 
strong-coupling region. As (3 increases, the action factor e~^ SQCU becomes relevant and 
will reduce the randomness of x^i 3 ) to some extent. Near the continuum limit (3 — > oo, 
however, the factor sin^(s) tends to approach in the MA gauge as shown in Fig.2(b), 
and hence such a constraint on x M (s) from the action is largely reduced, and the strong 
randomness of Xfi{ s ) is expected to hold there. Moreover, the reduction of the absolute 
value |c M (s) | = | sin # M (s) | itself further reduces the contribution of the off-diagonal element 
|c M (s)| in the MA gauge. 

Now, we examine the randomness of using the lattice QCD simulation. It should 
be noted that the residual U(l)3 gauge degrees of freedom should be fixed to extract 
itself, because x^( s ) is the U(l)3 gauge variant. To this end, we add the U(l)3 lattice 
Landau gauge JJB], which maximizes 



using the residual U(l) 3 gauge transformation in the MA gauge. In the U(l) 3 Landau 
gauge, there remains no local symmetry and the lattice variable becomes mostly contin- 
uous and approaches to the continuum field under the constraint of the MA gauge fixing. 
For the test of the randomness of x M (s), we calculate the probability distributions of x M (s) 
and the correlation between x^( s ) and X^( s + m the MA gauge with the U(l)3-Landau 
gauge. If Xfj,( s ) is a random angle variable, there is no bias on the distribution of x^( s ) 
and no correlation is observed between x^( s ) and X^( s + £)■ We show in Fig.5(a) the 
probability distributions P(Xn) and P(9u) at (3 — 2.4. Unlike P(du), P{x») is flat distri- 
bution without any structure in the whole region of j3, which is necessary condition of the 
random angle variable. We show in Fig.5(b) the probability distribution P(Ax) of the 
correlation 



which is the difference between two neighboring angle variables, at (3=0, 1.0, 2.4, 3.0. 
In the strong-coupling limit (3 = 0, Xti.{ s ) is a completely random variable, and there 



R — 22 -Retr 



(48) 



(49) 
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is no correlation between neighboring X\x- 111 the strong-coupling region as (3 < 1.0, 
almost no correlation is observed between neighboring Xni which suggests the strong 
randomness of Xfi( s )- As a remarkable feature, the correlation between neighboring x M 
seems weak even in the weak-coupling region as (3 > 2.4, where the action factor e _/3SQCD 
becomes dominant and remaining variables and 0^(s) behave as continuous variables 
with small difference between their neighbors as A9^ ~ and ~ 0. Such a weak 
correlation of neighboring x M would be originated from the reduction of the accompanying 
factor sin^ M (s) in the MA gauge. Moreover, in the weak-coupling region, the smallness of 
sin 9fj,(s) makes c M (s) more irrelevant in the MA gauge, which permits some approximation 
on Xfj.( s )- Thus, the random- variable approximation for x^( s ) would provide a good 
approximation in the whole region of (3 in the MA gauge. To conclude, the origin of abelian 
dominance for confinement in the MA gauge is stemming from the strong randomness of 
the off-diagonal angle variable Xni s ) anc ^ the strong reduction of the off-diagonal amplitude 
| sin 0^ (s) | as the result of the MA gauge fixing. 

VI. SUMMARY AND CONCLUDING REMARKS 

In the 't Hooft abelian gauge, QCD is reduced into an abelian gauge theory, and the 
color-magnetic monopole appears as the topological object in the constrained nonabelian 
gauge manifold corresponding to the nontrivial homotopy group n2(S77(iV c )/{7(l) iVc_1 ) = 
Z^~ x . Hence, if off-diagonal gluons can be neglected and the monopole is condensed, the 
QCD vacuum in the abelian gauge is described as the abelian dual superconductor and 
the confinement mechanism is understood as the dual Meissner effect. 

In relation with the dual Higgs picture for the confinement mechanism in QCD, we have 
studied the mathematical features of the abelian gauge fixing, the local gluon properties in 
the maximally abelian (MA) gauge, and the origin of abelian dominance for confinement 
in a semi-analytical manner with the help of the lattice QCD. 

First, we have studied the residual symmetry in the abelian gauge, with paying atten- 
tion to the global Weyl symmetry, which can remain as the relic of SU(iV c ). The global 
Weyl symmetry provides the ambiguity on the electric and magnetic charges, and persis- 
tently remains in the MA gauge. Considering the abelian gauge fixing in terms of the coset 
space of the fixed gauge symmetry, we have derived the criterion on the SU(A c )-gauge in- 
variance of the operator in the abelian gauge: if the operator defined in the abelian gauge 
is invariant under the residual gauge transformation, it is also SU(A c )-gauge invariant. 

Second, in the continuum SU(iV c ) QCD, we have expressed the MA gauge fixing us- 
ing the SU(A c )-covariant derivative operator. The local MA-gauge fixing condition and 
the composite Higgs field $[A M (s)] to be diagonalized are naturally derived from this 
expression of the MA gauge. 
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Third, we have examined the abelian projection rate i?Abeb the overlapping factor 
between SU(2) and abelian link variables, in the lattice formalism. In the MA gauge, we 
have found the high abelian projection rate as (i?Abei)MA — 1 for the whole region of f3, 
which means microscopic abelian dominance on the link variable. Using the statistical 
consideration, we have analytically estimated the lower bound of the abelian projection 
rate in the MA gauge as (-RAbei)MA ^ 0.84, which seems consistent with the lattice result 

(-RAbel)MA > 0.88. 

Finally, we have studied abelian dominance in terms of off-diagonal gluons in the 
Wilson loop in the MA gauge. In the SU(2) link variable, the off-diagonal angle variable 
is not constrained by the MA-gauge fixing condition at all, and tends to take random 
values besides the residual gauge degrees of freedom. By approximating the off-diagonal 
angle variable as a random variable, we have proved that the contribution of off-diagonal 
gluons to the Wilson loop, W Q s = (W / S u(2))/(W / Abei)MA, obeys the perimeter law in the 
MA gauge, which is numerically confirmed using the lattice QCD Monte Carlo simulation. 
Thus, we have showed exact abelian dominance for the string tension as 0su(2) = CAbei, 
although the finite size effect of the Wilson loop in the actual lattice QCD simulation 
leads to small deviation as <7su(2) > °Abei- 

In conclusion, we have found that the origin of abelian dominance for confinement in 
the MA gauge is stemming from the strong randomness of the off-diagonal angle variable 
Xp( s ) an d the strong reduction of the off-diagonal amplitude | sin# M (s)|, and these two 
remarkable features on the local variables Xfi( s ) an d O^s) are peculiar to the MA gauge 
fixing. 
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FIGURE CAPTIONS 



Fig.l: The gauge transformation property of $ and gauge function Q G G/H. The 
abelian gauge fixing is realized by Q G G/H so as to diagonalize $. (a) After the 
gauge transformation by v g G G, the operator <3> 9 is diagonalized by the gauge func- 
tion Q 9 = h[g]Qg^ G G/H. (b) The gauge transformation property of the operator O n 
defined in the abelian gauge. If O n is iY-invariant, Q is found to be invariant under the 
whole gauge transformation of G. 

Fig. 2: (a) The probability distribution P(i?Abei) of the abelian projection rate i?Abei at 
(3 = 2.4 on the 16 4 lattice from 40 gauge configurations. The solid curve denotes P(i?Abei) 
in the MA gauge, and the dashed line denotes P(Rxbe\) without gauge fixing, (b) The 
average of the abelian projection rate (i?Abei) in the MA gauge as the function of (3. For 
comparison, we plot also (Rxhei) without gauge fixing. 

Fig. 3: The off-diagonal gluon contribution on the Wilson loop, W oS = , as 

the function of the perimeter length L = 2(1 + J) in the MA gauge on 16 4 lattice with 
(3 = 2.4. The thick line denotes the theoretical estimation in Eq.(|40|) with the microscopic 
input (ln{cos6' M (s)})MA — —0.082 at f3 — 2.4. The data of the Wilson loop with I = 1 or 
J = 1 are distinguished by the circle. 

Fig. 4: The inter-quark potential V(r) as the function of the inter-quark distance r. 
The lattice data are obtained from the Wilson loop in the MA gauge on 16 4 lattice with 
(3 = 2.4 and T = 7. The square, the circle and the rhombus denote the full SU(2), the 
abelian and the off-diagonal contribution of the static potential, respectively. The thin 
line denotes the theoretical estimation in Eq. (|43|) . Here, the lattice spacing a is deter- 
mined so as to produce a = 0.89 GeV/fm. Due to the artificial finite-size effect of the 
Wilson loop, the off-diagonal contribution V oS gets a slight linear part. 

Fig. 5: (a) The probability distributions P(Xfi) (solid line) and P(0V) (dash-dotted 
curve) in the MA gauge with the U(l) 3 Landau gauge at (3 — 2.4 on the 16 4 lattice 
from 40 gauge configurations, (b) The probability distribution P(Ax) of the correlation 
Ax = mod vr (|x At (s) — Xpl(s + u)\) in the same gauge at (3 = (thin line), 1.0 (dotted 
curve), 2.4 (solid curve), 3.0 (dashed curve). 
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